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We consider the solution of a generalized E ® e Jahn- Teller Hamiltonian in the context of quasi- 
exactly solvable spectral problems. This Hamiltonian is expressed in terms of the generators of 
the osp(2, 2) Lie algebra. Analytical expressions are obtained for eigenstates and eigenvalues. The 
solutions lead to a number of earlier results discussed in the literature. However, our approach 
renders a new understanding of "exact isolated" solutions. 
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The Jahn- Teller (JT) distortion problem is an old one, dating back over sixty years 0. In 1958 Longuet Higgins 
et al. 0] carried out the first diagonalization of the Hamiltonian of the E ® e JT system within the adiabatic 
approximation. Yet, even today, new contributions to this problem are being made Q. They appear, however, not 
\ to have been fully exploited in the analysis of the JT problem. The E ® e JT problem is a system with a doubly 
■ degenerate electronic state and doubly degenerate JT active vibrational state. The JT effect describes the interaction 
of degenerate electronic states through non-totally symmetric, usually non-degenerate nuclear modes. This effect 
plays an important role in explaining the structure and dynamics of solids and molecules in degenerate electronic 
states. 

Studies of the JT effect led Judd to discover a class of exact isolated solutions of the model . To determine the 
relations between the parameters of the model, one can obtain an analytic form of two eigenvectors of the Hamiltonian 
corresponding to the specific energy. The complete description of these solutions has been given by Reik et al[f|. They 
ij observed that isolated solutions can be obtained by using Neumann series of expansions of the eigenvectors in the 
(Jj ■ Bargmann-Fock space described by the boson operators. The same problem has been treated by Klenner in polar 
coordinates Q. The possibility of the existence of some exact eigenvalues was previously proposed by Thorson and 
Moffitt . Their results have been generalized by Loorits || . 

As a related topic, the concept of quasi-exactly solvable (QES) systems discovered 0,^3 m tnc 1980s, has received 
much attention in recent years, both from the viewpoint of physical applications and their inner mathematical beauty. 
It turns out that in quantum mechanics there exist such some systems; a part of their spectrum can be computed 
using algebraic methods. The classification of the 2x2 matrix differential equations in one real variable possessing 
a polynomial solution has been carried out The relevant algebraic structure is the graded algebra osp(2,2). 

The E ® e JT Hamiltonian can be expressed as two coupled differential equations. The 2-dimensional QES matrix 
Schrodinger equation, with one real and one Grassmann variable have also been investigated|l2l Il5j. More recently, 
we have shown that the E ® (3 JT Hamiltonian is quasi-exactly solvable In this paper, we present a quasi-exact 
solution of the E (g) e JT Hamiltonian. 

This paper is organized as follows. In §2, following Reik's approach, we express the E ® e Hamiltonian in the 
Bargman-Fock space to obtain a second-order differential equation. In §3, this differential equation is expressed as 
linear and bilinear combinations of the generators of osp(2, 2) algebra. The condition for quasi-exact solvability of 
the E (8 e JT Hamiltonian is then derived. In §4, we discuss and compare our results with those studied given in the 
literature in different contexts. 



THE E ® e JAHN-TELLER HAMILTONIAN 



The well-known form of the E <g> e JT Hamiltonian describing a two- level fermionic subsystem coupled to two boson 
modes has obtained by Reik[5j is given by 

H = a+ ai + a+a 2 + 1 + + 2/i)cr° + 2/t[(ai + a+)cr+ + (a+ + a 2 )cr~], (1) 
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where \ + 2[i is the level separation, and n is the coupling strength. The Pauli matrices cr^ ,0 are given by 
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and they satisfy the commutation relations 

[*+,*-] =o°, [a°,a^} = T^. 
The annihilation and creation operators, a% and af, satisfy the usual commutation relations, 

[a+, a+] = [ai, aj] = 0, [ai,af] = <% . 
Following the analysis of Reik0], we can express the Hamiltonian in the form 

where J represents the angular momentum of the system and is given by 

J = a+ai — a^a 2 + -cr . 
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Note that J commutes with h, and the eigenvalue problem of the angular momentum part can be easily solved. It 
reads 



J |V) 



3 + 1 



(./ = <>. i.2 . . . ) 



with the eigenfunctions 



V>,- + i = (atYMafat) |0) |T) + (a+y +1 M44) |0> U) 



(7) 



(8) 



where |0) is the vacuum state for both bosons. Here ||) and ||) are the eigenstates of the operators a , and <j>i and 
cf>2 are arbitrary functions of afa^ ■ Because the operators h and J commute, the eigenfunctions |(SJ| are also the 
eigenfunctions of the Hamiltonian JQ). Therefore, we can write the eigenvalue equation 



and the equivalent Schrodinger equation 



whose eigenvalue e and E are related by 



H\i>) jH =E\^) jH 
3 

E = 2e + j + -. 



(9) 
(10) 

(11) 



From the above consideration, it is found that the solution of the Schrodinger equation can be reduced to the 
solution of H10JI . The Hamiltonian h can be expressed in the Bargman-Fock space by using the realizations of the 
bosonic operators 



i = 1,2. 



(12) 



In this formulation, the Hamiltonian h consists of two independent sets of first-order linear differential equations. 
Substituting (JHJ and H12|l into (|10[1 and defining ^ = zi.22 we can obtain the following two linear differential equations 
satisfied by the functions 4>i and 4>2- 



[f^-(e-M)]^i + «K^ + (S + i + l)]02 = 0, 
«[^ + l]0i + [f^-(c + Ai)]^a = 0. 



(13a) 
(13b) 
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These coupled differential equations represent the Schrodinger equation of the E ® e JT system in Bargmann's 
Hilbert space. In what follows, we introduce a new variable, 

i = n 2 (l + x), (14) 

and redefine the functions 

>h(0 = -~V>2(x), <l>i(£) = <pi(x)+<p2(x), (15) 
to obtain the following systems of differential equations: 

[-x-f- + (e + k 2 - fi)<pi(x) + \-x4- + 2e + 2k 2 + 1 + j + n 2 x]tp 2 (x) = 0, (16a) 
ax ax 

[4- + K 2 ]ipi{x) + \-x— + (e + k 2 + lAtpifx) = 0. (16b) 
ax ax 

It will be shown in the next section these equations can be solved in the framework of quasi-exact solvability. 

QES OF THE HAMILTONIAN 

The general procedure to solve a differential equation quasi-exactly is to express it in terms of the generators of a 
given Lie algebra having a finite dimensional invariant subspace and use algebraic operations. In a recent paper, [l6| 
we employed the sl 2 (R) Lie algebra in application to the E ® (3 system. In the present problem, the corresponding 
Lie algebra is a subalgebra of the osp(2, 2) Lie algebra. Let us introduce the following operators of the osp(2, 2) 
algebra (ljj: 

d 2 n , d _ , 

,7+ = X, J— = X— — — — Zk— hC (7 . 

ax z ax 

Jo = x— fc + -er~(7 + . J = k + —(T~a + , (17a) 

ax 2 2 

Qi = c~ ', Q2 = er~ — , Qx = —a + x- 1 — \-2ka + , Q 2 — xa + . (17b) 
ax ax 

Here J; (i = 0,+,—) are the bosonic generators and Qi,2 and Qi,2 are the fermionic generators. These generators 
form the osp(2,2) algebra. The bosonic generators form shiR) algebra, 

[J+, J-] = -2J , [Jo, J±] = ±J±, (18) 

and the fermionic generastors satisfy the anticommutation relations 

{Qi,Qj} = {Qi-.Qj} = 0, i = l,2. (19) 

The number operators of the system are given by 

N l = {Q 1 ,Q 1 }, N 2 = {Q 2 ,Q 2 }, (20) 
and the action of the number operators on the fermionic generator is defined as 

[JVi.Qi] = [N 1 ,Q 1 ]=0,[Ni,Q 2 }=Q 2 ,[N 1 ,Q 2 ]=-Q 2 , 

[N 2 ,Q 2 ] = [N 2 , Q a ] = 0, [N 2 ,Qi] = Q x , [N 2 ,Q X } = -Q x . (21) 

The matrix differential operators act on two-component spinors with components 

P NM = ({x°,x 1 ,...x N ,), x°a + ,x 1 a+,...x M a + ). (22) 

The linear and bilinear combinations of the operators given in (|17b|l Lie algebra is QES, when the space is defined 
as P n +i,n 15] . Let us consider the following linear combinations of the fermionic generators: 

L = 2/iiVi + (1 + 2fi)N 2 + k 2 {Q 1 + Q 2 ) + Q 2 + Q x . (23) 
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Then the eigenvalue problem can be expressed as 

L<p{x) = Xtp(x), (24) 

where 



tp{x) 



(25) 



is a two-component spinor. 

The algebraic structure of osp(2, 2) has been studied previously 0,^1- It is obvious that the fermionic generators 
of osp(2,2) with number operators Ni and N 2 form a subalgebra of osp(2,2). Inserting the differential realization 
(|17b(l into (|23[) and comparing the equations (|24|l with (|16a[) and (|16b() . we can show that the equations are identical 
when the following hold: 

A = i(l+j + 2/i + 2fc(l + 4^)), (26a) 

e = k - | - l - - k 2 . (26b) 

The restricted value of e are the same as those of well-known results 4J . When the generators act on the space 
Pn+i,n w e can obtain the two recurrence relations 

1 j 

nv n + (k - n + fi - - - -)u n + n 2 v n+l = 0, (27a) 
I j 

(2k - n)ui n + (k-n-fi-- - -)v n+ i + k 2 w„+i = 0, (27b) 

where uj n and v n are degree n polynomials in x and constitute a basis for the representation space of the generators 
of the algebra. They are related to the eigenfunctions as 

(P n+ i,P n ) = (ipi,(p2) = (v n+ i,uj n ). (28) 

It is necessary that the determinant of these sets be equal to zero, giving the compatibility conditions that establish 
the locations of the Juddian points on the energy baselines. The recurrence relation implies that the wavefunction is 
itself the generating function of the energy polynomials. The eigenvalues are then given simply by the roots of such 
polynomials. If Kj is a root of the polynomial (i> n +i, w„), the wave function is truncated at n = 2k and belongs to 
the spectrum of the Hamiltonian L. 



RESULTS 



The recurrence relations in (|27a|l and (|27b|l . depending on the choice of the parameters \x and j, describe a number 
of physical systems 0. These recurrence relations can be given in the matrix form: 



k + fi-^ 
2k 






K 
1 










1+2 
2 



K 

k + n - 2±i 
2k - 1 




k - /x - 
2 



Vl 
W2 



= 0. 



The matrix Hamiltonian is identical with the results obtained by Loorits [s| when we replace j by — j — 1 and set 
2/i = G. If /j, = and j is replaced by — j — 1, then our results are identical to Judd's results The physical systems 
described by the recurrence relations in (|27a|l and (|27b|) after replacing j by — j — 1, are summarized as follows. 

a) [1 = and j = 0: displaced harmonic oscillator. 

b) fi = and j is half integer: linear E ® e JT system. 

c) fi = and j is integer: linear T$ <S> t 2 and T% ® (s + t 2 ) systems. 

d) /j, ^ and j = 0: dimer. 

e) \x 7^ and j is half integer: linear E ® e system in an external field. 
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Thus we see that the method described in this article shows that quasi-exact solutions of a class of JT problems 
can be obtained as the roots of the determinant of a matrix of order 2k. The hypothetical case a) corresponds to a 
displaced one-dimensional harmonic oscillator for which the energy equation holds for all k: 

E= (2fc+0 -k 2 . (29) 

In the cases b) and c), which correspond to three octahedral JT systems, have common features. The form of all the 
determinants is the same for all three systems, provided that in the Fg <g> (e + T2) system, the two modes e and Ti 
are equally coupled to the electronic state Tg. The case j — 0, as the limiting case in the E ® e JT system, indicates 
relationship between E eg e JT and the Rabi system. It is known that the Rabi system and the E®f3 JT system occur 
in the resonant excitation of double molecules and dimers. The E (g /3 JT system possesses s^-R) symmetries 16J, 
which is a subgroup of osp(2, 2). 

Finally, in the presence of an external field, the energies of the generalized E eg e JT system can be obtained by 
solving the recurrence relations l|27a|l and (|27b|) . without any restriction. These recurrence relations form polynomials 
in K, and the first few of them are given by, 

Px(k) = r), 

P 2 ( K ) = 8ry K 2 -(p+l)(^-l), (30) 
P 3 (rt) = 12877K 2 -8(7](3r ? (p+l)-4)-4(p+l)) K 2 + W ( /0 + 2)(7 ? 2 -4), 

for the values k = 0, 1/2 and 1 , respectively. The roots of these polynomials belong to the spectrum of the Hamiltonian 
In the last equation, we have defined j = — (77 + p + 2)/2 and /i = (77 — p)/4. Any exact solutions that may be 
found (even if they are quasi-exact) would be useful in testing and improving the approximate and numerical results. 



CONCLUSION 



It is well known that exact solutions of the generalized E®e have a direct practical importance. We have presented 
a quasi-exact solution of the generalized E eg e JT system. The present work gives a unified treatment of some 
earlier works. The method given here can be extended to other JT and multi-dimensional atomic systems. A further 
interesting implication of the method is the existence of the relation between the E®e JT system and the two-photon 
Rabi Hamiltonian. 

The success of our analysis leads us the connection between isolated exact and quasi-exact solutions. We have shown 
that the E (g e JT system becomes quasi-exactly solvable with an underlying osp(2, 2) symmetry. The method given 
in this paper constitutes the most general Lie algebraic analysis of the JT interaction for E <g> e systems. We conclude 
that it is possible to study isolated, exactly solvable JT and other optical systems in the framework of quasi-exact 
solvability. In this way, isolated exact solutions of the JT Hamiltonians are put in the correspondence with solutions 
of the QES Hamiltonians. Finally, we mention a few things about further studies of the JT Hamiltonians. Instead 
of the tedious transformations given in §1, the fermion-boson realization of the osp(2,2) algebra can be constructed 
and then the E (g e JT Hamiltonian can be expressed in terms of the generators of the corresponding algebra. This 
construction may lead a simple calculation. 
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